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We show that interacting bosons in a periodically-driven two dimensional (2D) optical lattice 
may effectively exhibit fermionic statistics. The phenomenon is similar to the celebrated Tonks- 
Girardeau regime in ID. The Floquet band of a driven lattice develops the moat shape, i.e. a 
minimum along a closed contour in the Brillouin zone. Such degeneracy of the kinetic energy 
favors fermionic quasiparticles. The statistical transmutation is achieved by the Chern-Simons flux 
attachment similar to the fractional quantum Hall case. We show that the velocity distribution of 
the released bosons is a sensitive probe of the fermionic nature of their stationary Floquet state. 


It is well establishecPH^ that a one dimensional (ID) 
Bose gas with a strong short range repulsion (or equiva¬ 
lently small density) - the Tonks-Girardeau gas, exhibits 
features of weakly interacting fermions. Gan a similar 
phenomenon take place and be observed in dimensions 
larger than one? One example of statistical transmuta¬ 
tion in 2D is provided by the composite bosoi^icture of 
the fractional quantum Hall effect (FQHE)PEI. In this 
case the kinetic energy is totally degenerate due to the 
Landau quantization and the ground state is solel y d e- 
termined by the interactions. It was recently showrP^HIll 
that it is actually enough to have the kinetic energy de¬ 
generate along a line in the 2D reciprocal space (so called 
“moat”) to achieve the statistical transmutatioiP. The 
moat-like dispersion was d iscuss ed earlier in context of 
Rashba spin-orbit couplin^^^^^^K^, and for certain lat¬ 
tices with more than one site per unit cell and next near¬ 
est neighbors hoppin^ ^^ l ^^ l ^^ i 

In this paper we show that the moat dispersion may 
be found even in simplest lattices (e.g. square) upon a 
suitable periodic driving. The lowest Bloch band adi- 
abatically evolves into a Floquet band, which exhibits 
an approximately flat minimum along a closed contour 
encircling the F-point. Gonsequently the single particle 
density of states (DOS) diverges as inverse square root 
of energy at the bottom of the moat. Such ID behav¬ 
ior of 2D DOS motivates the analogy with the Tonks- 
Girardeau case. The specific mechanism of the statistical 
transmutation in 2D, however, is very different from its 
ID sibling. While the latter is achieved by simple sign 
inversion of the wave function at coinciding coordinates 
of particles, the former requires flux attachment, simi¬ 
lar to FQHE. As a result, the emerging 2D fermions are 
subject to the time reversal symmetry breaking effective 
magnetic field, and leading to a peculiar Landau spec¬ 
trum. We show that both the fermionic nature of the 
state and the effective magnetic field may be detected 
through the velocity distribution of the released gas. 

The Hamiltonian of interacting Bose gas with the moat 



c) 



FIG. 1: (Color online) a) Lowest Floquet band exhibiting ap¬ 
proximately flat moat, b) Frag ment of a square optical lattice 
with two sites per unit celt^ that gives raise to an approxi¬ 
mately flat moat in the lowest Floquet band upon resonant 
driving, c) Left panel shows the band structure of the un¬ 
driven optical lattice. The driving frequency uj is of order of 
the gap causing resonant coupling between the lowest and the 
first excited bands at k = 0. Middle panel: enlargement of 
the lowest and shifted first excited bands. Right panel: Blow 
up of the driven lowest Floquet band exhibiting a double-well 
feature in kx + ky direction indicating appearance of the moat 
in the Brillouin zone. 


dispersion relation can be written as 

^ i E N (|k| - ^o)Tr + 


( 1 ) 


where M is the mass of bosons, ko is the radius of the de¬ 
generate minima of the dispersion, g is the dimensionless 
contact interaction strength and ff>r are correspond¬ 
ingly boson creation and annihilation operators. 

To engineer a system of spinless bosons effectively de- 
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scribed by the Hamiltonian we adopt the idea of Flo- 
quet optical lattice shaking developed in Ref. |23| The 
authors employed harmonically shifted ID optical lat¬ 
tice at a near-resonance frequency corresponding to the 
lowest-band to the first-excited-band transition. As a re¬ 
sult they have created a ID dispersion relation with dou¬ 
ble minimum. The 2D generalization of this strategy is 
discussed below. A straightforward idea of using simple 
square lattice with the potential ^ sin^(/cLx) +sin^(/cL^) 
does not work because of the separable (in x, y) nature 
of the potential. Indeed the corresponding bands are la¬ 
beled by the two-integers €n^,ny with the lowest band 
being eo,o- The low energy excited bands eo,i are ei^o are 
highly anisotropic and thus do not lead to a flat moat. 
The next excited band is approximately isotropic, 
but shaking the lattice either in x or ^ directions does 
not produce matrix elements between and 60 , 0 ? due 
to orthogonality of the separable wave functions. 

The simplest way to engineer an approximately flat 
moat is to create a square optical lattice with two sites 
per unit cell. Such a lattice can be constructed by fusing 
a laser setup with laser beams of wavelength Xl = 27r jkL^ 
resulting in a standing wave intensity pattern forming a 
regular square lattice, and then adding two additional 
identical laser beams directed along x^y and x — y diag¬ 
onals (X scheme). If the lasers directed along x, and 
x — y directions are calibrated to have exactly the same 
phases while the laser along x -\- y diagonal has a phase 
shift amounting to tt, such setup results in a potential 


U{x,y) = Uo sin^ (kLx)+ 8111^ (kLU) 


This representation is analogous to the Bloch represen¬ 
tation of states in time independent periodic potentials. 
Here En{k.) G (0,co’) is the Floquet energy of n’th Bloch 
band, defined modulus of multiples of uj (corresponding 
to the energy quanta produced by driving ). The Bloch 
momentum is given by k = 27rm/I/, m = (ma,,m^), 
mx,y = with L being the lattice size. Here the 

function is periodic in x and y directions with 

the period and in time with the period T = 2t:Iuj. To 
find the Floquet spectrum we expand the periodic coun¬ 
terpart, (/9^^k(r, t), of the Bloch eigenfunction in Fourier 
series over momenta , ‘^Epv. Q^j^d energies ^ with 
integer s^mx^rriy G (— 00 , 00 ) and numerically diagonal¬ 
ize the Hamiltonian in this basis, assuming it is a large 
finite dimensional matrix. 

To understand the origin of the moat dispersion qual¬ 
itatively consider the two lowest bands, eo(k) and ei(k), 
of the potential 0. Their Bloch functions (/:?o,k(i*) and 
^i,k(i*) are correspondingly even and odd functions with 
respect to interchange of A and B sites of the unit cell, 
i.e. X y ^ —X — y. Figure Ic depicts these two bands 
for Uq = llF’i?, where Er = k\/2M is the recoil en¬ 
ergy. The resonant shaking of the phases with frequency 
UJ ~ ei(0) — eo(0) ~ 2.32Eji shifts the upper band as 
ei(k) — cj to touch the lower one at k = 0 and induces 
the off diagonal matrix element 14: between these two 
bands. Due to the parity of the Bloch functions, the lat¬ 
ter originates from the last term in Eq. (0 and is given by 
y(k) = ^Ji(4fcLA) / (fripo,ki^)sm[2kL{x + y)]<fi^]i{r), 
where Ji is the Bessel function. The effective two-band 
Hamiltonian 


+ sin^[/cL(x - y)] + cos^[/cl(x -h y)] 


( 2 ) 


This potential realizes a square optical lattice of the 
depth Uq^ having two sites per unit cell (A and B), 
with vectors /c^ei = (f, f), kLe 2 = (f:~^) and 
kL^s = (—^, 1^), depicted in Fig. 1 . One can har¬ 
monically drive this lattice by applying time dependent 
phase shifts to the lasers, resulting into the transforma¬ 
tion X ^ X — A coscjt, y ^ y — A coscjt, where A is the 
shaking amplitude. Below we show that, if the frequency 
UJ is nearly resonant with two lowest bands, an almost 
flat moat appears in the lowest Floquet band. 

The Lagrangian describing such a time-dependent 
problem is given by 


C = 




Hit) = -^ + U[x-xit),y-y{t)], (3) 


where '0(r, t) is a single particle wave function. Period¬ 
icity of lattice with = 'k jki,^ implies that the solution 
'0n,k(r,t) of the Shrodinger equation with Hamiltonian 
^(t), can be uniquely described by the space-time pe¬ 
riodic single particle eigenstate (/:?n,k(i*, t) of the Floquet 
operator idt — H(t) as ' 0 n,k(i*, t) = t). 


Hk 


eo(k) l/(k) \ 

C(k) ei(k)-w ) ■ 


( 4 ) 


yields the moat shape for the band, which is adiabati- 
cally connected to the original lowest Bloch band eo(k). 
Figure 0 shows dependence of the characteristic ra¬ 
dius of the moat as a function of Uq at two values of 
/cLA/ 7 r = 0.01 and 0.02. The moat appears at finite val¬ 
ues of UqIEr and exists until about Uq/Er < 14. Figure 
Eb shows the the ratio of the difference of maximal and 
minimal energies along the moat, Emax — ^min, and the 
hight of spectrum in the F point, Er — Emin, illustrating 
the flatness of the moat as function of Uq. 

An exactly flat Floquet moat band may be achieved 
by rotating of the laser^^ in ID optical lattices with the 
same resonant frequency as specified in Ref. m- 

Assuming an ideal moat, the interacting bosons are 
described by the Hamiltonian Q. Its single particle 
density of states diverges as similar to the ID 

case. This immediately suggests that the condensate 
can not be stable at any finite temperature. A more 
rigorous consideration. Ref. [25], starts from assuming 
a condensate in a state ko along the moat and derives 
the anisotropic spectrum of Bogoliubov excitations Ep = 

[e^ + 2(gn/M)ep]^^^, where ep = (|p+ko|-fco)^/2M. It 
is easy to see then that the fraction of thermally excited 
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FIG. 2: (Color online) Left panel: characteristic radius ko 
plotted as a function of dimensionless UqIEr at two values of 
— 0.01 (circles) and 0.02 (squares). Right panel: the 
flatness measure F of the moat versus UqIEr is shown. For 
definition see the main text. 
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quasiparticles diverges as a power law (Ref. [25] deals 
with 3D, where divergence is only logarithmic) with the 
system size, indicating instability of the condensate at 
any finite temperature. 

Moreover, there are compelling reasons to believe (see 
also Ref. [26|) that the groundstate of the repulsively in¬ 
teracting bosons with the moat dispersion does not break 
U{1) symmetry but breaks the time reversal symmetry. 
Thus the state is not represented by the condensate. It 
was argued in Refs. [T4UT6] that the better approximation 
to the actual ground state is the variational wave func¬ 
tion, given by a Landau level fully filled with fermions, 
which are transformed to bosons with the Chern-Simons 
(CS) flux attachment: 

$(ri, ■ ■ - riv) = • • -rjv). (5) 


Here ^^^(ri, • • • r^v) is a fully antisymmetric state of N 
fermions. By this reason the wave function (§ does 
not cost any short-range interaction energy. The CS 
phase = ( 2 ;^ — Zj)/\zi — Zj\, with complex 

Zj = Xj -h iyj, is antisymmetric with respect to the ex¬ 
change of any two coordinates, restoring the bosonic na¬ 
ture of ^(ri, • • • ttv). The CS factor costs kinetic energy, 
since it effectively modifies the momentum operator in 
Eq. 0 as iVr ^ iVr -A(r), where the vector potential 
is given by 


A„(ri) = Sap E 


(r* - rj)/3 


(6) 




and Cay is an antisymmetric tensor. The CS magnetic 
field, originating from A is given by B(Yj) = curl A(rj) = 
(5(ri — Vj) = 27rn(rj), where n(r) is the density 
operator. 

By analogy with the fractional quantum Hall effect 
we first adopt a mean field approximation, which re¬ 
places a set of the flux lines with a uniform mag¬ 
netic field B{r) B = 27rn, where n is the av¬ 

erage density. In this approximation is a state 
of non-interacting fermions with the dispersion relation 
(|k| — ko)‘^l2M placed in a uniform magnetic field B. 

The energies of corresponding Landau levels are given by 
r 1 2 


El = {kl/2M) 


{kl/2M) 


{I + 1 / 2 ) - 1 


where uor = 


FIG. 3: (Golor online) Excitation spectrum in units of Er is 
plotted versus njno. The per particle energy, EgsjER^ repre¬ 
sented by the fully filled lowest LL, is shown by the thick grey 
(red) curve. The thin black lines represent excited energies. 
Dashed line represents the chemical potential of a condensate 
Mconi/Ei^ corresponding to the interaction parameter g ^ 1. 
For such strong interactions the transition from composite 
fermion state to zero temperature BEG takes place at densi¬ 
ties n/no ~ 1. Upon lowering the transition shifts towards 
smaller n/no, as schematically shown in the inset. 


B/M = 2'KnjM and and I = 0,1,2,— The corre¬ 
sponding spectrum, shown in Eig. consists of set of 
non-monotonic functions of the field (i.e. density), which 
reach zero energy at 


ni 


47r(/-l-l/2)‘ 


(7) 


At this particular set of densities the fermionic wave func¬ 
tion is given by the fully occupied (indeed there is exactly 
one flux quanta per particle) l-th Landau level: 




—= det 
VA! 'rnd 



(8) 


where Xm(i*) = (27r(/+ m)!/!)“^/^(6’*‘)^+’^(a'*‘)^ ^ 

is a state with the angular momentum m = —I -/ + A 

at the Landau level L Here a’*' = = 

^(| — 2^) are the ladder operators. Eor intermediate 
densities < n < ni the ground state is a mixture 
of two fermion liquids with densities ni and The 

variational state 0,0 breaks time-reversal symmetry 
due to the effective Chern-Simons magnetic field, how¬ 
ever it conserves the U{1). Though within the mean- 
field approximation the state has zero energy, its actual 
energy is given by the expectation value of the Hamilto¬ 
nian operator 0 (or rather only its kinetic energy part). 
The corresponding calculations, given in the supplemen¬ 
tary material, show that the energy per particle scales as 

E{n) ^ ^ (^) This should be compared 

with either the naive estimate for the condensate state 
E ^ gn/M^ or the result of Ref. [26] E ^ In any 
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FIG. 4: (Color online) Dimensionless velocity distribution 
p{v) • vi^ where mvi = 7rni/ko, of an expanding gas of com¬ 
posite fermions plotted vs dimensionless velocity v/vi at fixed 
values of density n = / = 0, 3, 5 and zero temperature. In¬ 

set: Dimensionless velocity distribution of condensed bosons 
with density n = no. The temperatures are marked. At low 
temperatures the distribution shows a sharp peak at zero ve¬ 
locity, which is absent in the distribution of the composite 
fermions. 


event, the composite fermion variational state ( [^, Q is 
seen to be advantageous at small enough densit}^^ 

An important experimentally relevant measure of the 
composite fermion state of bosons discussed above is the 
velocity distribution of an expanding gas^hich can be 
observed in the time of flight experiments^. The group 
velocity of an expanding gas is defined by the derivative 
of the kinetic energy as v = ~ 

pectation value of this operator in the proposed state of 
composite fermions (§, (§ is obtained numerically and 
depicted in Fig. The result demonstrates striking dif¬ 
ference with the velocity distribution of condensed bosons 


shown in the inset of Fig. While at high tempera¬ 
tures distribution functions of condensate and of com¬ 
posite fermion state are similar, the qualitative differ¬ 
ence at T ^ 0 is caused by the fermionic nature of 
the latter. If for condensed bosons the distribution is 
sharply peaked at = 0, indicating condensation into a 
state with zero velocity, for composite fermions it is rem¬ 
iniscent to the Fermi-Dirac distribution exhibiting weak, 
plateau-like behavior at finite v at very low temperatures 
and small densities. Importantly, at low temperatures, 
there is no sharp peak at v = 0. The plateau vs peak 
difference can be regarded as the indication of the pro¬ 
posed statistical transmutation. In the field-theoretical 
language this difference can be traced back to the pres¬ 
ence of the effective Chern-Simons magnetic field and 
to the fact that effectively fermions find themseives in a 
state corresponding to the fuiiy occupied iowest Landau 
ievei. 

To conciude, we note that the abiiity to controi and 
probe statisticai transmutation in quantum many body 
systems is one of the most fundamentai chaiienges in 
contemporary physics. In this ietter we propose an ex- 
perimentai scheme to (i) engineer a resonantiy driven 
bosonic system exhibiting a moat-band and the phe¬ 
nomenon of transmutation of statistics, and (ii) probe 
our prediction for the veiocity distribution in time of 
flight experiments at low densities. The proposed state 
for bosons in a moat band is energetically more efflcient 
at low densities than any other known candidate for the 
ground state. It realizes a Floquet topological phase of 
bosons, joining the family of other topological structures 
achieved in non-equilibrium including the Floquet topo¬ 
logical insulators^^H^ and superfluid^. 

This work was supported by the PFC-JQI (T.S.), US- 
ARO, Australian Research Council, and Simons Foun¬ 
dation (V.G.), and DOE contract DE-FG02-08ER46482 
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i>(ri,---rjv) = H I + + F(ri,---riv) 


(SI) 


i<j 


where z = x-\-iy is the complex representation of the 2D vector r = (x, y), and I numerates the fully filled Landau level 
(LL). The wave function 4^i?(ri, • • is a fermion Slater determinant of spinless noninteracting fermions, residing 
on the l-th LL: 


• •-Tiv) = 




XmiZj) 


(S2) 


Here Xm (z) is the single particle wave function corresponding to a state with angular momentum m. It is given by 


Az) = 


{-iyVT\ 


^B\/27r2^(/ + m)\ 


Nr 


iO 


(m) 


I2-, 


[2^1 J’ 


m = -1. 


(S3) 


where L^^\x) is the adjoint Laguerre polynomial and Ib = 1/V27rn is the magnetic length. 

Non degenerate speetrum: As a first step we investigate the energy in the absence of the moat, ko = 0. It is given 
by the expectation value of the kinetic term (k^) (due to the antisymmetric fermionic multiplier of the wave function, 
the expectation value of the contact interaction potential is zero). In this expression, the Chern-Simons factor in (SI) 
creates a U{1) gauge vector potential, ^^(r) = One observes that, only one, two, and three particle states 

contribute to this expectation value. There are two types of terms contributing into (k^) = {k^)diag + 0^‘^)non-diag- 
contributions coming from diagonal in m, r variables states and non-diagonal ones. The general analytical expression 
for diagonal terms reads 


(k^)dms = /:Y / dr Xm(r){(k-.4(r))+?fdjag(r)}xm(r) 


(S4) 


where hereafter we suppress index I in most intermediate expressions and the gauge field is given by 

N-l N-l 


A^{Y)=e^^ Y (m,r'|H^(r-r')|m,r') = 


giv 


r',m=—Z 


Y /dr'x:„(r')-j^3/7yXm(r')- 


(S5) 


m=—l 
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The function 1-Ldiag{^) is given by 

'Hdiag{^) = T] /^J'i|Xmi(ri)pA 2 (r-ri) 

mi 

- E E{(”^i’ ri|(m2,r2|A(r -ri)A(r -r2)|m2,ri)|mi,r2) 


(S6) 


mi,m 2 ri,r 2 


+ ({mi,ri,m2,r2}|A(ri-r)A(ri-r2) + A(r2-ri)A(r2-r) |{mi, ri, m2, ri 


where the bra (ket) states (m,r| (|m,r)) represent the single particle function Xm{'^) (Xm (^))5 while the two-particle 
state |{mi, ri, m 2 , r 2 }) is given by a Slater determinant of two wave functions Xmi(i*i) and Xm 2 {'^ 2 ) with the same 
LL index 1. 

It is important to notice that in the thermodynamic limit the vector potential A^(r) gives raise to a constant 
Chern-Simons magnetic field 


A^(r)|jv^cx 3 ^ ^ B — 27m. 


(S7) 


Although the expression Eq. (S6) does not exhaust all contributions, the contribution of non-diagonal terms vanishes 
in the thermodynamic limit. Thus we conclude that the first term in Eq. (S4) exactly reproduces the mean-field 
result. 

The expression for the remaining non-diagonal matrix elements of (k^) can be cast as follows 

X^)non-diag =-^ E r| I {k, A(r - Ti)} + A2(r - Ti)|mi, r)|m, Ti) 


m,r,mi ,ri 


2 

N 


E E |^(m,r|(mi,ri|(m2,r2| |^A(r - ri)A(r - r2) + A(ri - r)A(ri - r2) 

r,ri,r 2 m,mi,m 2 

+ A(r2 - ri)A(r2 - r) |mi, r)|m2, ri)|m, r2)j, 


(S8) 


Where we used a notation {k. A} = kA + Ak. Due to the orthogonality of the single particle wave functions Xm (i*), 
most of this these non-diagonal terms are vanishing as A ^ 00 . Up to the terms vanishing in the thermodynamic 
limit A ^ oc, the non-diagonal terms yield 

{k'^)non-diag = -^ E (m, r| (mi, Fi | A^ (r - Fi) |toi , r) |to, Fi) 

= -^E/^"lx™WI^/<^rCTA^(F-Fi) 


. I 


Li 


r - ri 


2ll 


(S9) 


where we have used the completeness relation of the wave functions Xm 


EUHi’i)x(;i'Hr 2 ) = 




|ri -r 2 | 

21 % 


(SIO) 


We see that, using the completeness relation, one can cast the non-diagonal term (S9) in the form similar to the 
diagonal one, as in Eq. (S4). Joining Eq. ([S^ with 1-Ldiag{^) one obtains in the thermodynamic limit 




(Sll) 


where 




_ k-^i r 

A(r ——e L‘f 


|r - ri 


I2-. 


2 /| 


- / dYidY2 


27.2 


{2n)H% 


B Lf 


ri -F 2 


21 % 


(f-fi)2 + 27| 
A(f-Fi)A(f-F2) - 


(f-fi)2 + 27| 


ai 

2l%- 


(S12) 
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FIG. SI: (Color online) Dots represent the coefficient ai plotted vs the number I of the filled Landau level. The full line is the 
large I asymptote which is consistent with ai ~ 21og(7r/). 


Here we took into account the fact that last two terms in (S 6 ) give vanishing contribution due to the antisymmetric 

nature of the Slater determinant. We also have added and subtracted 7 - , 0/2 terms in both, first and second 

integrals which explicitly shows convergence at large distances in e ach o f them. Easy to see that the sum of the two 
added terms is identically equal to zero. The coefficient ai in Eq. (S 12 ) is position independent as its r-dependence 
may be absorbed into the shift of integration variables ri and r 2 . The coefficients ai can be calculated analytically 
by using the explicit form of Laguerre polynomials 




(S13) 


Substituting Eq. (S13) into Eq. (S 12 ), one is able to integrate the sum term by term. This results in the following 
series representation for the coefficient ai 


a/ — 1 + 


E 


/C1,/C2=0,/Ci+/C27^0 


(-i)''i+''^n 2 (fci + fc 2 )! 
{I - ki)\{l - k2)\h\2k2p 


H, 


/ci -|-/C2 — 1 5 


(S14) 


where = ^j=i j are the Harmonic numbers. Coefficient ai is depicted vs I in Eig.|^ With an excellent prescision 
ai is fit with the function 

ai ^ 21og(7r/). 

Energy in the presenee of a moat-band: In contrast to the case with non-degenerate dispersion, a moat-band 
spectrum studied in this paper exhibits a singularity. The \k\ singularity in the average kinetic energy (e(k)) = 
((|k| — ko)‘^) is however treatable. It can be dealt with by representing as 


= k^ 


f 


dse 




(S15) 


Using such a representation, one can straightforwardly exte nd th e expression (S4) to the case with a moat. One can 
see that the k dependence of an integrand in the rhs of Eq. (S15) is free from singularities. Thus, using the procedure 
discussed in the previous section one obtains 


1 


/' 


12 


((|k| - kof) = 77 E / V + ^o(r) - ko Xm{r) 


(S16) 


Assuming that the contribution to the energy coming from 'Ho{r) is mu ch sm aller than k^ (this assumption however 
needs to be checked self-consistently) and expanding the square root in (S16) over 'Ho(r)/fcQ, one obtains 

^o(r) 


((|k| 






drXmiO y(k-^(r)) - ko 


if (k-^(r))' 


+ ■ 


Xm(r)- 


(S17) 
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At densities ni = k^/ [47r(/ + 1/2)], when the mean-field yields zero energy measured from the bottom of the moat, one 

can separate the mean-field result from the fluctuation corrections. With the help of the identity ^(k — ^(r))^ |m, r) = 
/co|m,r) (which is the statement of the zero-energy within the mean-field) and Eq. (S12) the leading fluctuation 
correction is found to be 


E = ^m-ko)^) = 


1 


2MN 

1 

8Ml% 


?/**”<■■'(t) 


2ko) ~ 2 \Mk^ 


Xm(r) 

log^(4n/fco). 


(S18) 


Eq. (S18) constitutes the main result of the Sup plementary material. The validity of self-consistent assumption 
allowing to neglect higher order terms in Eq. ( S17) is g uaranteed when the scale ^ (^/^o) ^ smaller than one. 

This implies that the fluctuation induced energy (|S18|) is smaller than ^ n/M, which is satisfied precisely in the low 


density limit when the composite fermion state is more efficient than a condensate. 

If the number of particles is finite, the fluctuation effects give raise to Wdependent finite-size corrections to the 
energy. These corrections are vanishing in the thermodynami c lim it N ^ oc. Our analysis shows that the latter 
however reaches slowly. There are corrections to the energy (S18) originating from {^‘^)non-diag that scale as ^ 
n/{MN). There is also possibility of logarithmical corrections suggesting that in order to reach the thermodynamic 
limit numerically one has to take exponentially large systems, N > exp{/co/n}. 
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